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Abstract. - We generalize the invariant imbedding theory of the wave propagation and 
derive new invariant imbedding equations for the propagation of arbitrary number of coupled 
waves of any kind in arbitrarily-inhomogeneous stratified media, where the wave equations are 
effectively one-dimensional. By doing this, we transform the original boundary value problem 
of coupled second-order differential equations to an initial value problem of coupled first-order 
differential equations, which makes the numerical solution of the coupled wave equations much 
easier. Using the invariant imbedding equations, we are able to calculate the matrix reflection 
and transmission coefficients and the wave amplitudes inside the inhomogeneous media exactly 
and efficiently. We establish the validity and the usefulness of our results by applying them 
to the propagation of circularly-polarized electromagnetic waves in one-dimensional photonic 
crystals made of isotropic chiral media. We find that there are three kinds of bandgaps in these 
structures and clarify the nature of these bandgaps by exact calculations. 



Introduction. - The phenomena of the coupUng of two or more wave modes in inhomoge- 
neous media and mode conversion between them are ubiquitous in various branches of science, 
including plasma physics, optics, condensed matter physics and electrical engineering [1-6]. 
In this Letter, we develop a generalization of the powerful invariant imbedding method [6-13] 
to the case of several coupled waves in stratified media. Starting from a very general wave 
equation of a matrix form, we derive a new version of the invariant imbedding equations for 
calculating the reflection and transmission coefhcients and the field amplitudes. By doing 
this, we transform the original boundary value problem of coupled second-order differential 
equations to an initial value problem of coupled first-order differential equations. This makes 
the numerical solution of the coupled wave equations much easier. Furthermore, our equations 
have a great advantage that there is no singular coefficient even in the cases where the material 
parameters change discontinuously at the boundaries and inside the inhomogeneous medium. 
We check the validity and the usefulness of our invariant imbedding equations by applying 
them to the propagation of electromagnetic waves in stratified chiral media. By calculating 

(*) E-mail: khkimOaj ou . ac . kr 
© EDP Sciences 



2 



EUROPHYSICS LETTERS 



the matrix reflection and transmission coefficients exactly, we clarify the nature of the three 
different photonic bandgaps that can exist in photonic crystals made of chiral media. 

Theory. - We consider a system of N coupled waves propagating in a stratified medium, 
where all parameters may depend on only one spatial coordinate. We take this coordinate as 
the z axis and assume the inhomogeneous medium of thickness L lies in < z < L. We also 
assume that all N waves propagate in the xz plane. The x component of the wave vector, 
is a constant and the dependence on x of all wave functions can be taken as being through a 
factor e*''^. In a large class of interesting problems, the wave equation of N coupled waves in 
the present situation has the form 

S - fz^~'^'^fz + [^(^)^'-^(^) - - 0, (1) 

where ij) = ('01, • ■ • , V'Af)"^ is an A^-component vector wave function and £ and A4 are N x N 
matrix functions that depend on z in an arbitrary manner inside the inhomogeneous medium. 
We assume that the waves are incident from the vacuum region where z > L and transmitted 
to another vacuum region where z < 0. / is a unit matrix and K is a diagonal matrix such that 
Kij = kiSij, where ki is the magnitude of the vacuum wave vector for the i-th wave. £ and 
M are unit matrices in the vacuum region. The nonsingular functions £{z) and M{z), which 
specify the material properties of the medium and/or the external conditions, can change 
discontinuously at the boundaries and at discrete z values inside the medium. By assigning 
£{z) and A4{z) suitably, eq. is able to describe many different kinds of waves in a large 
number of stratified media. 

There are numerous examples where the effective wave equations have precisely the same 
form as eq. Later in this Letter, we will apply our theory to the propagation of electromag- 
netic waves of two different polarizations in layered chiral media, where "0 is a two-component 
vector and £, A4 and K are 2x2 matrices. Another interesting example is the propagation 
of the probe and phase-conjugate waves in layered phase-conjugating media [14, 15]. A wide 
variety of mode conversion phenomena observed in space and laboratory plasmas can also be 
studied using eq. 1^ [1,2,13,16]. 

Following Gryanik and Klyatskin [6] , we generalize eq. slightly, by replacing the vector 
wave function by an iV x iV matrix wave function the j-th column vector (^ij, • • • , ^nj)"'" 
of which represents the wave function when the incident wave consists only of the j-th wave. 
We are interested in the N x N reflection and transmission coefficient matrices r — r{L) and 
t = t{L). Let us introduce a matrix 

where T and T are the time-ordering and anti-time-ordering operators respectively. P is a 
diagonal matrix satisfying Pij ~ PiSij and pi is the negative z component of the vacuum wave 
vector for the i-th wave. It is straightforward to prove that g{z, z') satisfies the differential 
equations 

d d 

-T^giz, z') = i sgn(z - z') £{z)Pg{z, z'), -^giz, z') = ~i sgn(z - z') g{z, z')£{z')P. (3) 
oz oz' 

Using eqs. (|2l and the wave equation ^ is transformed to an integral equation 
*(z,L)=g(z,L) 

- ^ ^ dz'g{z, z') [£{z')P - PM{z') - q^p-'M{z') + q^p-^£-\z')\ L), (4) 
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where we consider as a function of both z and L. We take a partial derivative of this 
equation with respect to L and obtain 

^^^=i^S{z,L)a{L) + <i>{z,L), (5) 

where 

a{L) = £{L)P - i [£{L)P - PM{L) - q^p-^M{L) + q^p-^£-^{L)] L), (6) 

and $(z, L) satisfies an equation similar to eq. I@J except that there is no source term (that 
is, g{z, L)). This implies <i>(z, L) — Q and then we have 

'-%^=^n^,L)aiL). (7) 
Taking now the derivative of 'i'{L, L) with respect to L, we obtain 



d^{L,L) _ d'^iz,L) 



dL dz 



a*(z,L) 



dL 

=L 



= i£{L)P[r{L)- I]+i^{L,L)a{L). (8) 

=L 



Since '^{L,L) = / + r{L), we easily find the [N x N matrix) invariant imbedding equation 
satisfied by r{L): 

fif 

— = ^[r{L)£{L)P + £{L)Pr{L)] 

-\[r{L) + /] [£{L)P - PM{L) - q^p-^M{L) + q^p-^£-\L)\ [r{L) + /]. (9) 

Similarly by setting 2 = in cq. Q, we find the invariant imbedding equation for t{L) 
(=*(0,L)): 

^ = ^t{L)£{L)P 

-U{L) [£{L)P - PM{L) - q^p-^M{L) + q^p-^£-\L)] [r{L) + I]. (10) 

These invariant imbedding equations are supplemented with the initial conditions, r(0) = 
and t{0) = I. For given values of P and q and for arbitrary matrix functions £{L) and A4{L), 
we solve the coupled nonlinear ordinary differential equations (O and H1U|) numerically using 
the initial conditions, and obtain the reflection and transmission coefficient matrices r and t 
as functions of L. The invariant imbedding method can also be used in calculating the field 
amplitude ^'(z) inside the inhomogeneous medium. Rewriting eq. (UJ, we get 

-^*(z, [£{l)P - PM{1) - q^p-^M{l) + q^p-^£-^{l)\ [r{l) + I]. (11) 

For a given z (0 < z < L), the field amplitude ^{z, L) is obtained by integrating this equation 
from I — z to I — L using the initial condition ^(z, z) — I + r(z). 
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Application. - Eqs. ©, H1U|) and which have never been derived before to the 

best of our knowledge, wiU be the starting point in our future analysis of a variety of wave 
coupling and mode conversion phenomena. In the rest of this Letter, we establish the validity 
and the utility of our invariant imbedding equations by applying them to the problem of the 
electromagnetic wave propagation in stratified chiral media. 

Isotropic chiral media are those where the appropriate constitutive relations are given by 



D = eE + i7H, B = - i^yF,. 



(12) 



The parameters e, fi and 7 are the dielectric permittivity, the magnetic permeability and the 
chiral index respectively [17-19] . Some researchers use alternative constitutive relations [20,21] 



D = eE + i^B, H = B/^ + «^E. 
The two relations give identical results if the parameters are identified by 



(13) 



(14) 



We will use eq. (|12|) from now on. In recent years, there have been a large number of theoretical 
[17, 19-24] and experimental [18, 25-28] studies on the wave propagation in various kinds of 
chiral media. 

From the Maxwell's equations and the constitutive relations, we are able to derive the 
wave equations satisfied by the electric field in inhomogeneous chiral media: 



1, 



/uV X ( X E ) (e/i - 7^) + ^ 



7V X E - 



/xV X ( -^E 



(15) 



In the uniform case, right- and left-circularly-polarized waves are eigenniodcs of this equation 
with the effective refractive indices ^^fjl+j and ^JT\x—^ respectively. In inhomogeneous media, 
these two modes are no longer eigenmodes and are coupled to each other. The equation 
satisfied by the magnetic field H is similar except that the roles of e and [i are reversed. 
In media stratified in the z direction, e, [i and 7 are functions of z only. For plane waves 
propagating in the xz-plane, the x dependence of all field components is contained in the 
factor e**^. In this situation, we can eliminate -E^,, i?^, Hx and from eq. (|15|l and obtain 
two coupled wave equations satisfied by Ey = Ey(z) and Hy = Hy{z), which turn out to have 
precisely the same form as eq. with 



Ey 



K 



k 
k 



£ 



fi «7 
-17 e 



M 



e 



(16) 



where k = lo /c. 

We have used eqs. ®, (IIUII and (|16|) in calculating the reflection and transmission coeffi- 
cients in various situations. In all cases where exact solutions by other methods are available, 
our theory gives the same results. In our notation, rii(r2i) is the reflection coefficient when 
the incident wave is s-polarized and the reflected wave is s(p) -polarized. Similarly, r22(?'i2) is 
the reflection coefficient when the incident wave is p-polarized and the reflected wave is p{s)- 
polarized. Similar definitions are applied to the transmission coefficients. By a suitable linear 
combination of these coefficients, we are able to obtain a new set of the reflection and trans- 
mission coefficients r^j and , where i and j are either -f or — [19]. For instance, r_|_+(r y.) 

represents the reflection coefficient when the incident wave is right-circularly-polarized and 
the reflected wave is right(left)-circularly-polarized. The reflectances and transmittances are 



defined by Ri 



fijl'^ and Tij 
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As an example, we consider a uniform chiral layer of finite thickness with the parameters 
e, ^ and 7, placed between uniform achiral media of infinite thicknesses. In this case, the 
electromagnetic wave equations can be solved analytically, following the methods used in 
elementary quantum mechanics. Lekner has presented an exact analytical solution of this 
problem [19](^). In defining the reflection and transmission coefficients, rss, ''sp, ''ps, ''pp, 
tss, tsp, tps and tpp, Lekner uses different conventions from ours. In order to compare his 
solution with ours, we need to identify rn = rss, '^12 = ^ps, r2i = —rsp, ^22 = — "^pp, tn = tss, 
ti2 = tps, t2i = tsp and t22 = tpp, if the magnetic permeability of the media outside the chiral 
layer is equal to 1. Taking these into account, we have verified analytically that Lekner's 
expressions for the reflection and transmission coefficients satisfy our invariant imbedding 
equations exactly. 

In fig. n we plot the refiectances R++, and and the transmittances T++ and 

T , when a wave is incident at 9 = 45° on a one-dimensional photonic crystal made of 

alternating chiral and dielectric layers of the same thicknesses A/2. It can be proved easily 

that i?i2 = i?2i and = R |-. The chiral layer has the parameter values of e = 4, = 1 

and 7 = 0.3 and the dielectric layer has e = 2, ^ = 1 and 7 = 0. The total number of periods 
is 50. The x component of the wave vector, q, is given hy q = Losm9/c and the z component 
of the vacuum wave vector matrix, P, is given by P = pi, where p — lucos9/c. Also plotted 
is the imaginary part of the Bloch wave number k for an infinitely large photonic crystal. 
This quantity was obtained using an exact analytical expression for the dispersion relation of 
infinitely large photonic crystals made of two different kinds of alternating chiral layers, which 
we have derived recently [29] . The frequency region where the imaginary part of k is nonzero 
corresponds to a photonic bandgap. 

We find an excellent agreement between the analytical result on the dispersion relation and 
the reflectance and transmittance spectra. In general, there are three kinds of bandgaps, two 
of which are so-called co-polarization bandgaps and one of which is called a cross-polarization 
bandgap. Unlike in previous studies of this phenomenon [23,24], our theory is free of any 
approximation and provides exact band structures. For large values of 7 and 9, these three 
bandgaps can be well-separated, as demonstrated in fig. ^ where we show the second group 

of bandgaps. The R-\ spectrum clearly displays a cross-polarization bandgap and the R++ 

and R spectra show co-polarization bandgaps. The transmittance spectra show that a 

right(left)-circularly-polarized wave, the frequency of which lies in the co-polarization bandgap 
of left{right)-circularly-polarized waves, is freely transmitted. 

It is straightforward to apply our method to more general situations where the parameters 
e, /i and 7 are arbitrary functions of z. For example, we can easily study the effects of defects 
and randomness on the wave propagation in chiral media using eqs. ©, ((11])) . ffTT^ and 
Our equations can also be applied to the cases where both e and fj, take negative values with 
no modification. A detailed study of this so-called negative refractive index medium [30], 
which is also chiral, is of great interest and will be presented elsewhere. We have also applied 
our method successfully to a number of other coupled wave problems, such as the phase- 
conjugate reflection of light from nonlinear phase-conjugating media, the light propagation in 
uniaxial and biaxial media and the mode conversion phenomena in both unmagnetized and 
magnetized plasmas [31]. In these studies, we solve the full wave equations exactly, without 
using common approximations such as the slowly varying envolope approximation and the 
WKB approximation. These results will be presented in a near future. 

(^)There are three typos in Lekner's solution. In the expressions of and in eq. {A2), cf + c+c— and 
+ C-\-C— have to be replaced by Cj — c+c_ and — c+c— respectively. The expression Z+Z^ appearing at 
the end of the equation for tps in eq. (A5) has to be replaced by Z'^Z—. 
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Fig. 1 - (a) Imaginary part of the Bloch wave number k for an infinitely large photonic crystal 
made of alternating layers of chiral and dielectric materials of the same thicknesses A/2. The chiral 
layer has the parameter values of e = 4, ^ = 1 and 7 = 0.3 and the dielectric layer has e = 2, 
= 1 and 7 = 0. The transverse component of the wave vector, g, is given hj q = ujsmO/c, 
where 6 = 45°. The z component of the vacuum wave vector matrix, P, is given by P = pi, 
where p = uicos9/c. The frequency region where Im k is nonzero corresponds to a bandgap. (b-f) 
Reflectance and transmittance spectra for a one-dimensional photonic crystal made of alternating 
layers of chiral and dielectric materials. The parameter values and the values of 9, p and q are the 
same as in (a) and the number of periods is 50. 
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